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ABSTRACT 



Context. The solar corona is heated to high temperatures of the order of 10^^. The coronal energy budget and specifically possible 
mechanisms of coronal heating (wave, DC-electric fields, ..) are poorly understood. This is particularly true as far as the formation of 
X-ray bright points (BPs) is concerned. 

Aims. Investigation of the energy budget with emphasis on the relative role and contribution of adiabatic compression versus current 
dissipation to the formation of coronal BPs. 

Methods. Three-dimensional resistive MHD simulation starts with the extrapolation of the observed magnetic field from SOHO/MDI 
magnetograms, which are associated with a BP observed on 19 December 2006 by Hinode. The initial radially non-uniform plasma 
density and temperature distribution is in accordance with an equilibrium model of chromosphere and corona. The plasma motion is 
included in the model as a source of energy for coronal heating. 

Results. Investigation of the energy conversion due to Lorentz force, pressure gradient force and Ohmic current dissipation for this 
bright point shows the minor eff'ect of Joule heating in comparison to the work done by pressure gradient force in increasing the 
thermal energy by adiabatic compression. Especially at the time when the temperature enhancement above the bright point starts to 
form, compressional eff'ects are quite dominant over the direct Joule heating. 

Conclusions. Choosing non-realistic high resistivity in compressible MHD models for simulation of solar corona can lead to un- 
physical consequences for the energy balance analysis, especially when local thermal energy enhancements are being considered. 

Key words. Sun: atmosphere - Sun: magnetic topology - Magnetohydrodynamics (MHD) - Methods: numerical - Sun: corona 



1. Introduction 

The mechanisms of coronal heating are not well understood. 
A particular object for studying heating processes are coronal 
bright points (further abbreviated BPs). Due to the increasing 
accuracy of observations our knowledge about BP has greatly 
advanced from the time of their discovery in soft X-ray images 
[Vaiana et al. | (|1970[). According to X-ray and EUV observa- 



tions the linear size of BPs is on average about 30-40 arcsec 
with, typically, an embe dded bright core of about 5-10 arcsec 
[Madjarska et al. (2003 ). The average lifetime of X-ray BPs is 
about 8 hours IGolub et alTl i P^T?) ) and 20 hours for EUV BPs 
[Zhang et alT ( |2001 ). For a long time it has been known that BPs 
are associated with small bipolar magnetic features in the pho- 
tosphere Krieger]et^7]([T97TJ^ About one 
third of BPs lie over emerging regions of magnetic flux, while 
the rest of them lie above moving magnetic features. This was a 
base for the "cancelling magnetic feature" (CMF) model Priest| 



et al. ( 1994). Lifetime and energy release of BPs are known to 



be closely related to the diflTerent phases of the motion of this 
photospheric magnetic feature |Brown et al7]p001| ). First theo- 
ries were mainly addressing t he topology of the magnet ic field 
below BPs by e.g., Parnell etlT] ([T953)); ^LongcopeT] ([1998]). 
Using higher resolution and cadence observations of BP's inten- 
sity and taking into account a more comprehensive patterns of 
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motion in particular in regions with highly divergent magnetic 
field, ( Brown et al. ||200T ) could associate diflTerent patterns of 
motion of the solar photospheric magnetic features to diflTerent 
stages of a BP evolution. The plasma motion in the regions of 
strong magnetic field was first included by Biichner (2004a,b) 
in their three-dimensional numerical resistive MHD model us- 
ing their 3D numerical simulation model, LINM0D3d. The lat- 
ter considers dissipation of currents generated by plasma mo- 
tion in photosphere on time scales longer than an Alfven time 
as a one of the heating processes in the solar corona [Parker, | 
(i l972| ). In their model they took int o account current dissipa- 



tion due to anomalous resistivity (Biichner & Elkina 2005 



200 6) that causes Joule heating. Since LINM0D3d considers 
the compressibility of the plasma, the resulting heating could 
be due also to compressional eflTects. Later o n two-dimensional 
MHD simulation studies were carried ou t byjvon Rekowski B. 
jet al. I ([2006a'b'); 'von Rekowski et al. (2008 ). These authors 
used an analytical initial equilibrium and imposed a magnetic 
flux footpoint motion to model coronal bright point heating as 
being due to canceling magnetic features. To obtain the desired 
heating rate they used an enhanced resistivity for which the val- 
ues were above the theoretically justifiable resistivity. This raises 
the general question of the energy budget and energy conversion 
in solar flux tubes. Even with low resistivity, current simulations 
are unable to resolve the diflTusion regions of reconnection and 
thus overestimate Joule heating. It is also unresolved how much 
heating is caused by pressure gradient forces. 
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To clarify this question we continued t h e work of Biichner 
P0Q4a|b|c|); |B iichner ' ( |2006l |200T); "Santos & Biichner* 



^ 2QQ7| )7 [Santos et al. (2008 ). These authors demonstrated the 
formation of localized current sheets in and above the transition 
region at the position of a EUV BPs as a result of photospheric 
plasma motion. This study is extending their results through a 
systematic study of the energy conversion and budget in mag- 
netic flux tubes. The investigation uses the 3D simulation model 
LINM0D3d to simulate the solar atmosphere in the region of an 
X-ray BP observed by the Hinode spacecraft on 19 December 
2006 between 22.17 UT and 22.22 UT. 

In section[2]we briefly review the main features of the numer- 
ical simulation model LINM0D3d. In section [3] we describe the 
specific simulation setup used in our study and section [4] pro- 
vides some simulation results for the chosen BP data. In sec- 
tion |5] we present results of energy budget analysis by investi- 
gating the role of diff'erent forces and in section [5| we summarize 
and discuss our results. 



2. Simulation model 

Our simulation model uses the approach of the LINM0D3d code 
( Biichner et al. |[20()4a|b|c ). This means that the initial magnetic 
field is obtained by extrapolating the observed photospheric line- 
of sight (LOS) magnetic fields. The initial plasma distribution is 
non-uniform containing a dense and cool chromosphere as well 
as the transition to a rarefied and hot corona. The photospheric 
driving is switched on by coupling the chromospheric plasma 
with a moving background neutral gas. Some details of our code 
have been given briefly in the following subsection. 



2.1. Equations 

In our study we solve the following set of MHD equations: 
dp 

dt 
dpu 

"dT 

dB 

~dt 

dp 



= -V • pu 
-V • puu - Vp -\- j X B - yp(u - Mo) 

= V x(uxB - T]j) 
-V -pu-iy- \)pV • M -h (r - \)rif 



(1) 
(2) 
(3) 
(4) 



where p and u are plasma density and velocity, B is the mag- 
netic field and P is the thermal pressure. A plasma-neutral gas 
coupling in photosphere and chromosphere is included through 
the collision term in the momentum equation, where Mq denotes 
the neutral gas velocity. The neutral gas serves as a frictional 
background to communication photospheric footpoint motion to 
the plasma and magnetic field through frictional interaction. It 
also leads to a reflection of coronal Alfven waves back to the 
corona from the transition region, so that the influence of coro- 
nal Alfven waves can be neglected at the photospheric boundary. 
In order to set the plasma in motion a number of incompressible 
flow eddies is used according to observed horizontal drifts in the 
photosphere V • uq = is imposed via the neutral gas, where Uq 
is dependent in x and y. It is constant along z and derived from a 
potential using uq = V x (f/e^), with 



U = uqo/ cosh 



x-y + cp 



I cosh I 

\ ^1 



y + ci 



(5) 



t/oo /^i , chosen in accordance with the observed plasma motion 
in the photosphere. In our simulation we approximated the ob- 
served motion by three vortices with amplitudes of the velocity 
Woo equal to 5.5, 5 and 2 hmj s, respectively. The values of co, 
Lo, ci and Li are 9, 6, 51 and 6 Mm for the first vortex 5, 6, 28 
and 6 Mm for the second and 19, 7, 38 and 7 Mm for the third 
vortex. The height-dependent collision frequency v is chosen to 
be sufficiently large only below the transition region. This way 
the plasma is forced to move, dragged by the neutral gas, in the 
model chromosphere but not above the transition region. This 
way the horizontal motion generates a Poynting flux into the 
corona. On the other hand the collision frequency is chosen in a 
way that coronal Alfven waves are properly reflected while wave 
perturbations in the chromosphere are heavily damped by the 
frictional interaction with the neutral background. Our choice 
of equations means that in this study we do not consider en- 
ergy losses due to radiation and heat conduction and we also 
excluded the action of the solar gravitation in this study. The 
system of equations is closed by Ohm's and Ampere's laws and 
the temperature is defined via the ideal gas law for a fully ionized 
plasma: 



E = -ux B -\-rjj 
VxB = fioj 
p = IukbT 



(6) 
(7) 
(8) 



The value of the resistivity ri is varied in accordance with 
three models described in subsection 2.3 The MHD equations 
are discretized by means of a second order weakly dissipative 
Leapfrog scheme. Due to stability reasons the induction equation 
is discretized using Dufort-Frankel scheme, , Potter , ( ,1973^ ). 

2.2. Simulation box and normalization 

The lower boundary of the simulation box is a horizontal square 
in the photosphere sized 46.4 x 46.4 Mm^. The simulation box 
extends 15.45 Mm toward the corona. A nonuniform grid in the 
z direction supplies the proper resolution of the transition layer, 
where the grid distance Az corresponds to 160 km, ^Biichner et| 
al. ( 2004a| ). This corresponds to 64 grid points in z direction, 
while in the x, y plane a 128 x 128 grid are used. We solve for 
dimensionless variables that are normalized to natural scales as 
listed in table. [T] Note that the maximum imposed velocity of 
the neutral gas is smaller than 5 km/s while the typical (nor- 
malizing) electron thermal velocity is Vthe = 1470 km/s and the 
Alfven speed is va = 50 km/s. Hence, one can be certain that the 
inserted neutral gas motion is gentle, sub-Alfvenic and sub-slow 
velocities. 



2.3. Resistivity models 

In order to verify the influence of difl'erent resistivity models on 
the BP plasma heating we solved the equations for the same ini- 
tial and boundary conditions but varying the resistivity model. 
The resistivity r] can be expressed via an efl'ective collision fre- 
quency yu as 7/ = - 



where ojpe is the electron plasma fre- 



Note that the contour lines of this function are streamlines 
of the flow. The magnitudes of velocity scale with uoq/Lq and 



quency (ojpe = ^Ine^ le^nie). In our model we always apply a 
constant physically justified background resistivity t/o which ex- 
ceeds exceeds the numerical resistivity. It is appropriate to chose 
for eff'ective col lision frequency of the background resistivity the 

I I 1 3 II 

Spitzer, ( 1962 ) value p = (ne'^LnAT'^ /ISne^m^ ). Based on 
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the typical plasma parameters of our model we chose for the 
collision-driven background resistivity 7/0 = 10""^ (in normalized 
units). In two models we switched on additional, anomalous, re- 
sistivity in places where either the current density of the current 
carrier velocity (uccv determined as the current density divided 
by the charge density) exceeds a physically justified thresholds 
of micro-instabilities. 

In the first resistivity model anomalous resistivity is switched on 
when the current carrier velocity (Uccv exceeds a critical velocity 
dRoussev et aTl [20021 [Blichner& Elkinal[2Q05] ): 



Table 1. Normalization values. 



^ = ^0 + 



0, 



if |w, 



ccv\ ^ l^crit 
crit 



(9) 



A natural choice for the threshold velocity is the electron thermal 
velocity Vthe^ in our for the normalizing quantities 1470 km/s or 
to 5.8.10"^ in normalized units. In the first resistivity model we 
chose 5.10"^ to follow the ideal evolution of the plasma as long 
as possible. The additional term for resistivity can be estimated 
e.g., fo r a nonlinear ion- acoustic instability ( Biichner & Elkina 
|l [2Qn6l )as 



neff 



(10) 



Here ojpi denotes plasma ion frequency (ojpi = ^jne^leQnii). For 
the typical parameters of our simulation this estimate would re- 
veal 7] = 2.5, i.e. a magnetic Reynolds number of less than 
unity. In this case many current sheets would immediately dif- 
fuse away. On the the other hand, since the plasma /3 is relatively 
large for our simulation parameters obliquely propagating waves 
would be present in the spectrum of the micro- turbulence. In this 
case the estimate of the eff'ective collision frequency has to take 



into account lower-hybrid waves ( Silin & Biichner 2005 ). For 



our normalizing values this results in rje/f = 0.03. 

In a second model calculation we considered a current den- 
sity dependent resistivity used before, e.g., by |Neukirch"et 



al. (1997), in which the resistivity increases even stronger 
(quadratic dependence) after the current density exceeds a criti- 
cal value jcrit- 



0. 



if\j\ < 7c, 



rieffif-^ - if\j\ > jcrit 



(11) 



The critical current density is related to critical velocity via 
jcrit = e He Merit- Here we will report the results of our simu- 
lations obtained according to the second model for which we 
chose a threshold as low as jcrit = 0.69 in order to discuss the 
consequences of an early switch on of additional, anomalous re- 
sistivity. For comparison we solved the problem also by assum- 
ing for a third model a constant enhanced uniform resistivity as 
usually done in global MHD simulations. 

Concerning the values of the chosen rje/f one should note 
that the width of the actual current sheets in which turbulence ef- 
fectively operates is of the order of the ion inertial scale ^ = — . 
This scale cannot be resolved in any realistic 3D MHD simula- 
tion of the solar corona. In order to introduce micro-turbulent 
anomalous resistivity the threshold velocity (- current density) 
has to up- scaled to the actual resolution of the simulation by a 
factor of 5.101 By the same reason the resistive electric field 
builds up in very (perhaps J/-) thin current sheets. To consider 



Variable 



Normalization value 



density 
lenght 

magnetic field 



Lo = 500 km 

5o = IG = IQ-"^ T 



pressure 
temperature 
Alfven velocity 
time 



Po= =4x10-3 //m2 



To = = 7.2 X 10^ K 

^ ^ riQ Kb 



VaO = 



To = Lo/ vao = 10 s 



= 50 km/s 



the correct values of the electric field on the much coarser MHD- 
simulation grid anomalous resistivity used in the simulation has 
also to be scaled up by the above scaling factor. This approach 
allows to consider the correct amount of Joule heating. 

2.4. Initial and boundary conditions 

We first carried out a potential field extrapolation to the Fourier 
decomposed normal field component of the magnetic field taken 
from the MDI-observation. The resulting 3D magnetic configu- 
ration is used as the initial condition of the simulation code. In 
the potential field approximation the normal field component is 
related to (B;,, By) through V • 5 = and V x 5 = 0. The initial 
density and temperature height profiles for the plasma is taken in 
accordance with the VAL model that assumes pressure being in 
a hydrostatic equilibrium. The simulation box has 6 boundaries: 
4 lateral, 1 top and 1 bottom boundaries. For the side boundaries 
a line symmetric boundary condition is used with the line sym- 
metry with respect to the centers of the sides of the simulation 
box. For the upper boundary the derivatives in the normal direc- 
tion are put to zero. At the lower boundary the normal velocity 
is set to be zero, while the tangential velocity is taken from the 
neutral motion. 



3. Simulation setup 

Our study is based on an X-ray BP observed by the XRT X-ray 
telescope on board of the Hinode spacecraft on 19 December 
2006. The corresponding X-ray image is shown in Fig. [T] For 
the initial magnetic field we used the observed line - of - sight 
(LOS) component of the photospheric magnetic field taken by 
the Michelson Doppler Interferometer MDI onboard the Soho 
spacecraft at 22:17 UT. For that sake data from a field of view 
with the horizontal size of 64 x 64 arcsec^ was chosen that prop- 
erly covers the magnetic features associated to this BP (insert in 
Fig.[T]). Note that we use the LOS component as the initial nor- 
mal field component at the lower boundary of our simulation 
box, the photosphere, since the BP observation was made close 
to the center of the solar disc. 

Fourier filtering was applied to the LOS component of the 
magnetic field. By taking into account only the first eight Fourier 
modes, details of magnetic field structure smaller that 6 Mm are 
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Fig. 1. X-ray Image taken from XRT/Hinode on 19 December 
2006 at 22:08 UT. Insert: the LOS component of the photo- 
spheric magnetic field in a 64 x 64 arcsec^ horizontal plane 
taken from MDI/Soho, where white (black) spots correspond to 
upward (downward) directed Line-of- sight components of the 
photospheric magnetic field. The BP and the related magnetic 
field feature are indicated in the images. 



neglected. The extension of structures arising from smaller scale 
magnetic features would not extend higher up into the corona, 
they are dissipated at an early stage of the evolution in the highly 
collisional chromospheric plasma. 




-40 -30 -?n -10 10 20 



Fig. 2. Potential magnetic field extrapolated from the filtered 
MDI magnetograms and used as initial configuration for the 
magnetic field in our simulation. The blue lines show the mag- 
netic field lines. The color code depicts the LOS component of 
the magnetic field. Note that axes here are in terms of grid points, 
64 in z direction and 128 x 128 in x,y plane. 

Fig. [2] shows a three-dimensional view of the magnetic field 
extrapolated from the photospheric boundary for the magnetic 



field observed at 22:17 UT on December 19, 2006. The blue 
lines show the magnetic field lines. The color code depicts the 
LOS component of the photospheric magnetic field. Magnetic 
fields directed upward from the photosphere are colored in red, 
downward directed in blue. 

With the chosen normalization length of Lq = 500 km, the 
box size in x and y direction correspond to 92.8 Lq and the z di- 
rection extend to 30.9 Lq. The photospheric plasma velocities 
are obt ained by applying the loca l-correlation-tracking(LCT) 
method [November & Simon 1 ( |1988| ) to the Fourier filtered LOS 
magnetic component of the photospheric magnetic field ob- 
served between 22:17 UT and 22:22 UT. The left panel of Fig.|3] 
shows the velocity pattern obtained by the LCT method. For the 
simulation we used incompressible velocity vortices to mimic 
the observed velocity pattern, as shown in the right panel of 
Fig. |3] Note that the interval chosen for the simulation starts a 
few hours after the time the BP first appeared in the X-ray im- 
ages and that the bright point continues to glow a few more hours 
afterwards. During the whole simulation time interval the rela- 
tive shear motion of the two main magnetic flux concentrations 
of opposite polarity is negligible. 







LCT time: 




22:17 to 22:22 UT 

























Fig. 3. Horizontal plasma velocities in the photosphere. Left 
panel: velocities obtained by applying LCT technique to MDI 
magnetograms between 22:17 and 22:22 UT. The right panel 
shows the vortices that are used to approximate this motion in 
the simulation. 



4. Simulation results 

The simulation results are first shown in a plane at x = 45.7 
(Fig.|2]), which crossed through the center of the two main mag- 
netic polarities. The vertical profile of the temperature is shown 
in Fig.|4]for t = (top panel), 80 (middle panel) and 160 s (bot- 
tom panel). In ^ = we have a height dependent temperature as 
defined by the initial condition. At ^ = 80 s the eff'ects of plasma 
compression and expansion, together with Joule heating, shape 
the temperature profile. An arc of hot plasma is formed above 
the two opposite magnetic polarities. The increase in tempera- 
ture in this layer is approximately 0.5 in normalized units, what 
corresponds to 36000 K. The region that is located just below it. 
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however, experiences some drop in temperature. At ^ = 160 s the 
arc of hot plasma leaves the simulation box and we are left with 
a corona in which the differences in temperature can reach one 
orders of magnitudes. 



J perpendicular B 



Temperature 



(x= 45.7) 




Fig. 4. Temperature distribution in the vertical plane at x = 45.7 
at the beginning (upper panel) and at t = 80 s and t = 160 s. 
Note the temperature in color bar in presented in terms of the 
normalization value. 

Fig. |5] shows the parallel and perpendicular components of 
the current with respect to the magnetic field direction at t = 80 
and t = 160s. It can be seen that enhanced current flows coincide 
well with the temperature increase. This would lead to an inter- 
pretation of the heating as being due to current dissipation only. 
However, as shown later, adiabatic heating can have an impor- 
tant contribution to temperature increase. 

5. Energy balance 

Let us now diagnose the difl'erent contributions to plasma heat- 



ing in the BP region. First, in subsection 5.1 we discuss the 



overall global heating. In subsection |5.2| the dependence on the 
resistivity model is presented. Finally, the flux-tube heating is 
analyzed in subsection [53 



5.1. Global effect of current dissipation and compression 

In order to understand the relative contribution of current dis- 
sipation and plasma compression to the coronal plasma heat- 
ing in the BP region it is appropriate to analyze the pressure 
changes rewriting the Eq. 4 in terms of a continuity equation for 
the Temperature evolution. This leaves two source terms on the 
right hand side of the equation: 

^ + V . = -(r- 1) r V . M + (r- 1) (12) 

Let us first analyze the first case, where an anomalous resistivity 
is used when the current carrier velocity exceeds a critical value. 




Fig. 5. Parallel and perpendicular components of electrical cur- 
rent at t = 80 s (upper two panels) and t = 160 s (lower two 
panels) in the same diagnostic plane as in Fig.4. Note the en- 
hancement in perpendicular current is located at the same place 
of the temperature maximum. 



Fig. [6] shows the resulting distribution of -(y - 1) T V • m in the 
vertical diagnostic plane. This way we have a proxy for tempera- 
ture changes associated to pressure compression and expansion. 

Adiabatic heating has an important role on the formation of 
the high temperature arc that propagates upward towards the top 
boundary. It is also due to expansion that temperature decreases 
below this hot arc. 

The values of second term in the right hand side of the Eq. 
11, (y - 1) rjf/ p, is shown in Fig. [t] in the plane x = 45.7 at 
t = 80s and t = 160s. By comparison with the compressional 
part the contribution of the Joule heating appears to be negligi- 
ble. For a better comparison the contribution of the two terms in 
the right hand side of the Eq. 1 1 in the temperature evaluation, 
the horizontal view is shown at the height of transition region in 
two diff'erent instance of time, t = 80 s in left and t = 160 s in 
the right panel of Fig. [8] In the following, we will analyze in 
some more detail to what degree compression and Joule heating 
contribute to the evolution of the temperature. For this sake and 
in order to study the role of the forces involved in the energy 
conversion process, we performed a volume integration of the 
time rates of change of kinetic, magnetic and thermal energies 
in the simulation box above the chosen Bright Point region. Our 
approach is similar to |Bim et al. | ( |2009] ) when they used energy 
transport equations to analyze the properties of energy conver- 
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Fig. 8. Temperature, first and second terms of Eq. 12 over density, in the left, middle and right panels, respectively. The result is 
shown at a horizontal plane in transition region at two instance of time, t = 80s in top and t = 160s in the bottom panels. 
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Fig. 6. Adiabatic cooHng/heating rate after 80 s and 160 s in the Fig. 7. Joule heating rate, (second term in the r.h.s of Eq. 12, 



plane x = 45.7, according to the first term in the r.h.s of Eq. 11 
over density, (-(y - l)TV • u <0). 



divided by density), after 80 s and 160 s in the plane x = 45.7. 



sions associated with a reconnection process. The contribution 
of diff'erent terms in the energy transport process can be studied 
from the following equations: 



dt 



= { pu^U'ds+ f(-U'Vp-\-U'jxB)d^v (13) 

2 Jsv Jv 



dSrr, 



dt 



dsth 
dt 



-r 



i-uB^ + (u • B)B - 7]j xB)'ds (14) 

+ I (-U • j X B - j]f)d\ 
Jv 



I pu • ds -\- I (u 



Vp + j]f)d'v (15) 



Where Skm, ^mag and Sth denote kinetic, pw^/2 , magnetic, 
B^/2fio, and thermal, P/(y - 1), energies, respectively. The vol- 
ume integrals (second term on the right-hand side) in these equa- 
tions represent the energy conversion from one form into the 
another. This energy conversion are explicitly written in terms 
of the work done by Lorentz force, pressure gradient force and 
Joule dissipation, (left panel of Fig. [9]). Note that the initial spike 
in the Lorentz force is in part caused by numerical discretization 
errors and in part by the onset of photospheric footpoint motion. 
The initial oscillations are damped substantially during, approx- 
imately, two Alfven times, followed by a state of an approximate 
force balance. This eff'ect was found to be smaller in a run where 
footpoint motion was not included. The initial perturbation has 
a minor eff'ect on the initial extrapolated magnetic field, it does 
not aff'ect the currents and Lorentz forces at a later times. 
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The surface integrals are also needed to obtain the energy 
rates, when they indicate the transport of each of the three form 
of energies. With the chosen boundary condition however, the 
values of these surface integrals are zero at the lower bound- 
ary. They compensate each other through the side boundaries of 
the simulation box as well. At the upper boundary however, one 
needs to consider the contribution of this surface integrals in the 
rate of energy transfer. This means E x B, Pu and pu^u for the 
transport of magnetic, kinetic and thermal energies, respectively. 
The values of these terms at the upper boundary are shown in 
Fig. 
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One can see that the contribution due to these terms is 
insignificant, so it would be a good approximation to consider 
only the volume integrals for the change in the energy rates. 
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Fig. 9. Rates of the contribution to the total energy changes, (bot- 
tom panel). Diff'erent contributions to the changes of the energy, 
(top panel). 



The changes in energy rates are shown in the right panel of 
Fig. |9] the forces responsible for these changes are depicted in 
the left panel of the Figure. As one can see by comparing the two 
panels the magnetic energy is transferred to kinetic energy al- 
most completely via the work done by the Lorentz force that ac- 
celerates the plasma. It is an intermediate step however, followed 
by the work done by pressure gradient force that converts kinetic 
energy into thermal energy. This decelerates the plasma motion 
until, finally, the Lorentz force is balanced. The direct transfor- 
mation of magnetic energy to thermal energy (Joule heating) is 
via Ohmic current dissipation, tjJ^. A comparison of the en- 
ergy conversions rates (see Fig. [9j right panel) however shows 



that Joule dissipation plays a minor role in the energy exchange 
process while the other contributions are orders of magnitudes 
larger. The minor role of Joule heating in comparison to adia- 
batic process in the increase of thermal energy was also found 
for the case of a solar flare by |Bim et al. ( |2009j , where they ex- 
plained the compressional heating in two almost simultaneously 
steps: acceleration by Lorentz force and deceleration by pressure 
gradients. 
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Fig. 10. First term in the right hand side of the Eqs.12-14, (sur- 
face integrals) is shown in the upper boundary of the simulation 
box after 160 s. 
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5.2. Influence of different resistivity models 

The previous calculation was based on an anomalous resistiv- 
ity model with the current carrier velocity as a critical value for 
a local switch-on of additional resistivity. In order to better un- 
derstand the influence of the resistivity we performed the sim- 
ulation also with two other resistivity models, one that uses a 
current density dependent resistivity and another with constant 
resistivity respectively. 

Fig. [TT] depicts the resulting energy conversion rates and the 
work done by the involved forces v • / x 5, v • VP and by 
J] for all the three resistivity models by using difl'erent line 
styles for the results obtained by using the difl'erent resistivity 
model. The results obtained for the three cases show that the 
resistivity model influences the dynamics of the system and the 
thermal energy rate mainly through the pressure gradient force. 
While magnetic and kinetic energy rates of change depend only 
weakly on the resistivity model, the rate of temperature change is 
significantly influenced. Nevertheless, independent on the used 
resistivity model the heating is due mainly to the work done by 
the pressure gradient force. At the same time the contribution 
of the Joule heating is about two orders of magnitude smaller 
(Note the scale of the plots in the top row.) We conclude that 
the adiabatic compression is the dominant eff'ect in increasing 
temperature in the BP region in all three cases. 



5.3. Flux tube heating 

In order to locate the heating eff'ect better it is appropriate to 
determine it for individual flux tubes, integrating along the mag- 
netic field lines instead of taking values averaged over the whole 
simulation box as reported in the previous sections. In this inte- 
gration one has to take into account the changing cross- section 
of flux tubes. This can be done by applying the concept of the 
difl'erential flux tube volume V = J B'^ds, where ds indicates 
the step size along the field line. This way the flux conservation 
in a flux tube (O = A x 5 = const.) is taken into account by 
the proportionally of the cross-section to B~^. Note that large 
flux tube volumes correspond to field line rising high into the 
corona or hitting regions of vanishing magnetic field. The energy 
is transported in accordance with the upward directed Poynting 
flux E X B, enhanced magnetic tension is carried away by wave 

propagation. 

For the quantities described in secti on |5.1| the resulting flux 
tube integrated values are shown in Fig. |12|in the horizontal ref- 
erence plane just above the transition region. The values reached 
indicate once more the negligible role of Joule heating by cur- 
rent dissipation for the thermal energy change in the bright point 
region compared to the dominant role of the pressure gradient 
force. Please note the difl'erent range of the plots in Figure 12 as 
indicated by the color bar. It also can be seen that the locations at 
which this force and also maximum rates of energy changes ap- 
pear coincide. Furthermore, the same pattern has formed in the 
integration result of v • VP, y - JxB and the rate of change of the 
difl'erent kinds of energy. This pattern can clearly be seen in the 
integration of total energy along the field lines (Fig. 13_), which 
is the sum of the kinetic, magnetic and thermal energies: 



Total energy 



pu^ -h —B^ -h -^}(Pv 
2 2yUo r - 1 



The left panel of Fig. 14 shows the result of this integration 



for temperature and flux tube volume. The coincidence of the 
temperature enhancement with the maxima obtained in the flux 
tube integrated energy change rates and forces shows that the 
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Fig. 13. Results of integration along the magnetic field lines us- 
ing the diff'erential flux tube volume concept for total energy(top, 
left panel) and magnetic(top, right panel), kinetic(bottom, left 
panel) and thermal energies(bottom, right panel), at t = 160s. 

heat is provided by the plasma compression due to the Lorentz 
force. 



Temperature 
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Fig. 14. Temperature and flux tube volume, integrated along the 
magnetic field lines, at t = 160s. 



Enhanced flux tube integrated values follow the same pat- 
tern as the BP. This indicates that the regions of enhanced tem- 
peratures correspond to the foot points of field lines leading to 
higher altitudes or to regions where the magnetic field vanishes. 
The plasma motion across these regions supplies the magnetic 
energy that is converted to thermal energy. 



6. Summary and discussion 

We have presented the results of heating processes in the region 
of an observed X-ray coronal bright point. In particular we have 
investigated the importance of the work done by adiabatic com- 
pression in comparison with Joule heating in the course of the 
dynamic evolution and heat production near the bright point. 

The simulation shows that an arc- shaped structure of en- 
hanced temperature forms that is 2-4 times hotter than the back- 
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Fig. 11. Energy change rates for three different resistivity models, (top panel). The bottom panel shows the work done by the 
Lorentz force, pressure gradient force and the Joule heating power, different lines correspond to anomalous current carrier dependent 
(dashed), anomalous current dependent (dotted) and constant (solid line) resistivity models. 
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Fig. 12. Integration along the field lines using the differential flux tube volume concept for the work done by Lorentz force, pressure 
gradient force and Joule heating, (top panel, from left to right), and the changes in rates of magnetic field change, the temperature 
and the kinetic energy, (bottom panel, from left to right) at t = 160s. 



ground plasma. This structure is located above the two main 
opposite photospheric magnetic flux concentration. It coincides 
with the location where the electrical current densities are max- 
imum. The structures of temperature and current density en- 
hancements, indeed, coincide. 

We further examined the contribution of the Lorentz force, 
pressure gradient force and Joule heating performing volume in- 
tegrals in the simulation box that determine the magnetic, kinetic 



and thermal energy change rates for three different resistivity 
models. We found that independent on the resistivity model mag- 
netic energy was transformed to kinetic energy through the work 
done by Lorentz force. Kinetic energy in turn is converted to 
thermal energy due to pressure gradients that balance the Lorentz 
force. 

A comparison of the effect of the three energy conversion 
through Y ' J X B,Y ' VP and 7/ show that adiabatic compres- 
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sion has an important role in temperature increase in the upper 
corona. This is not dependent on the resistivity model used in 
the simulation. 

For a better understanding of the heating processes we uti- 
lized the concept of differential flux tube integration of the dif- 
ferent contributions along the magnetic field lines. A quantitative 
comparison in the horizontal plane, from where the integration 
starts, shows that energy conversion rate, total energies and work 
done by Lorentz and pressure gradient forces are located in the 
same flux tubes, also temperature and flux tube volume are max- 
imum at the same place. 

We conclude that the conversion of magnetic energy to ki- 
netic energy via the work done by the Lorentz force and from 
kinetic to thermal energy due to the work done against the pres- 
sure gradient force determine the heating of this bright point. 
We could show that plasma compression dominates the heat- 
ing of the bright point. In contrast, the role of Joule dissipation 
appeared to be negligibly small. The temperature enhancement 
follows the same pattern. The fact that the pattern obtained by 
calculating flux volume integrals coincides with the one of tem- 
perature and energy change rates bring us to the conclusion that 
plasma motion at the footpoints of the flux tubes carries the en- 
ergy upward and makes the flux tubes rise to the higher corona. 
The magnetic energy is converted to thermal energy until the 
plasma compression is balanced by the Lorentz force. In the lo- 
cal, flux-tube oriented consideration we also could see that the 
role of the Joule heating in these energy conversion processes 
was negligible and the heating of plasma in the bright point re- 
gion is basically due to pressure gradient force. 

First, the fact that Joule heating is weak in the corona 
was not entirely unexpected but it is quantitatively con- 
firmed here. It is worth to remember that the necessary 
up-scaling of the resistivity and of the onset condition of 
micro-turbulent anomalous resistivity to the resolved by the 
MHD simulation grid scales does even overestimate the ac- 
tual Joule heating. As a result Joule heating cannot be con- 
sidered a viable process unless there is a convincing argu- 
ment that the dissipation regions are volume filling to a much 
larger extend than the already large one used in the present 
model. 

Second, the results demonstrate very clearly that compres- 
sion is an important processes in the energy budget. It is 
not clear in how far compression can contribute to the over- 
all coronal heating but it is certainly important for the local 
heating of BPs. 

Third, in this context the nature and the consequences of 
plasma compression are worth some consideration. In ideal 
MHD adiabatic compression is reversible. But the conse- 
quent fiux tube heating is, however, irreversible due to mag- 
netic reconnection and other mixing processes. Magnetic re- 
connection, in particular, changes flux tube identities (mag- 
netic connectivity) while fiux tube entropy conservation re- 
quires ideal MHD in addition to appropriate boundary con- 
ditions. Local adiabatic compression becomes irreversible 
also due to other plasma transport processes like heat con- 
duction and radiative cooling. These aspects will be sepa- 
rately investigated in a subsequent paper. Meanwhile the re- 
sults presented here clearly demonstrate that in the overall 
energy budget plasma compression (and expansion) can play 
an important role in the heating of the corona. 
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